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ABSTRACT 

The Serpens SMM 1 region was observed in the 6.9 mm continuum with an angular resolution of about 0'.'6. Two 
sources were found to have steep positive spectra suggesting emission from dust. The stronger one, SMM la, is the 
driving source of the bipolar jet known previously, and the mass of the dense molecular gas traced by the millimeter 
continuum is about 8 Mq. The newly found source, SMM lb, positionally coincides with the brightest mid-IR source in 
this region, which implies that SMM lb is yet another young stellar object. SMM lb seems to be less deeply embedded 
than SMM la. SMM 1 is probably a protobinary system with a projected separation of 500 AU. 

Subject headings: ISM: individual (Serpens SMM 1) — ISM: structure — stars: formation 



1. INTRODUCTION 

The Serpens dark cloud is a nearby star-forming region (see 
Harvey et al. 2007 and references therein). The cloud core con- 
tains a cluster of low-mass young stellar objects, and some of 
them are extremely young protostars (Testi & Sargent 1998; 
Harvey et al. 2007; Winston et al. 2007). The core also harbors 
several Herbig-Haro objects and molecular outflows (Davis et 
al. 1999; Hodapp 1999). The most luminous and most deeply 
embedded object among them is SMM 1 (Harvey et al. 1984; 
Casali et al. 1993; Enoch et al. 2007). SMM 1 is a Class source 
and associated with a bipolar jet/outflow (Rodriguez et al. 1989; 
Curiel et al. 1996; White et al. 1995; Larsson et al. 2000). 

Because of the extremely high extinction, the small scale 
structure of SMM 1 has been accessible only with radio in- 
terferometric observations. Images of centimeter continuum re- 
vealed a well-collimated bipolar jet showing a large proper mo- 
tion away from the central protostar, indicating a dynamical age 
of -60 yr (Rodriguez et al. 1989; Curiel et al. 1993). Hoger- 
heijde et al. (1999) imaged SMM 1 in the millimeter continuum 
and found a single compact source surrounded by an envelope 
of ~9 Mq. Hogerheijde et al. (1999) also found that the molec- 
ular outflow driven by SMM 1 has a complicated structure. 

One of the interesting issues in the study of star forma- 
tion is the multiplicity, because the majority of stars belong to 
multiple-star systems (Duquennoy & Mayor 1991). The mul- 
tiplicity of young pre-main-sequence stars is higher than that 
of main-sequence stars (Duchene et al. 2007), which implies 
that multiple systems form in the early stage of star forma- 
tion. In fact, bright protostars, such as IRAS 16293-2422 and 
NGC 1333 IRAS 4A, are usually found to be multiple sys- 
tems when imaged with a high angular resolution (Looney et 
al. 2000). Serpens SMM 1 has been considered as an interest- 
ing example of an isolated protostar that can be understood well 
because it is relatively bright, massive, and simple. However, 
there were some indications of complexity, such as the near-IR 
excess (Larsson et al. 2000). Therefore, it is necessary to inves- 
tigate the structure of SMM 1 with a high angular resolution and 
a high sensitivity. 

In this paper, we present the results of our observations of the 
Serpens SMM 1 region in the 6.9 mm continuum with the Very 
Large Array (VLA) of the National Radio Astronomy Observa- 
tory. We describe our radio continuum observations and archival 
data in Section 2. In Section 3 we report the results of the contin- 
uum imaging. In Section 4 we discuss the star-forming activities 
in the SMM 1 region. 



2. OBSERVATIONS 

The Serpens SMM 1 region was observed using VLA in 
the standard Q-band continuum mode (43.3 GHz or A = 6.9 
mm). Twenty-three antennas were used in the D-array config- 
uration on 2001 December 22, and twenty-six antennas were 
used in the C-array configuration on 2004 March 27. The phase- 
tracking center was a 2 ooo = 18 h 29 m 49 i ;65, S 200Q = 01 o 15'20'.'6. 
The phase was determined by observing the nearby extragalac- 
tic source 1851+005 (G33.50+0.20). 

In the D-array observations, the flux was calibrated by setting 
the flux density of the quasar 1331+305 (4C 30.26) to 1.18 Jy, 
which is the flux density measured within a day of our observa- 
tions according to the VLA Calibrator Flux Density Database. 
Comparison of the amplitude gave a flux density of 0.65 Jy for 
1851+005. In the C-array observations, the quasar 0410+769 
(4C 76.03) was observed for the purpose of flux calibration, but 
its flux around the time of our observations is not listed in the 
database. Since the flux of 0410+769 varies between 0.25 and 
0.55 Jy with a variability time scale of about a month, it could 
not be used for the calibration. The flux of 1851+065 is much 
more stable than 0410+769, and the flux scale of the C-array 
data was calibrated by setting the flux density of 1851+005 to 
0.65 Jy, as found in the D-array observations. The flux scale 
was checked by comparing the total flux densities of SMM 1 
from the two observing tracks. In this comparison, the D-array 
dataset was limited to make the shortest uv length as the same as 
that of the C-array dataset. The resulting flux densities of SMM 
1 agreed within 1%. Therefore, the flux scale of the C-array data 
is consistent with that of the D-array data. 

2.1. Archival Data 

In order to investigate the source structure and the contin- 
uum spectra, we analyzed several VLA datasets of the X-band 
continuum (8.5 GHz or A = 3.5 cm) and C-band continuum (4.9 
GHz or A = 6.2 cm) observations retrieved from the NRAO Data 
Archive System. The data were calibrated and imaged following 
a standard procedure. The NRAO observing project AC 504 in- 
cludes data from the A-array observations made on 1998 May 9, 
May 30, and June 1. The 3.5 cm map produced by us is consis- 
tent with that of Raga et al. (2000). However, the flux densities 
of the outflow sources are much weaker than what were seen 
eight years before (Curiel et al. 1993). Therefore, to measure 
the flux densities, we also analyzed the dataset of the NRAO ob- 
serving project AC 563, which includes data from the A-array 
observations made on 2000 October 27 (about 1 .2 yr before our 
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Fig. 1. — Uniform- weight maps of the 6.9 mm continuum {thick con- 
tours) toward the Serpens SMM 1 region from the observations in 2001 
December. The rms noise is 0.2 mjy beam -1 . The contour levels are 
1, 2, 4, and 8 times 0.6 mjy beam -1 . The 6.9 mm sources are labeled. 
(a) Map restored with the usual synthesized beam shown in the bottom 
right corner: FWHM = 2'.' 5 x l'/5 and P.A. = -45°. Thin contours: 
Map of the 6.2 cm continuum from the AC 504 dataset. The rms noise 
is 0.012 mjy beam -1 . The contour levels are 1, 2, 4, 8, and 16 times 
0.04 mjy beam" 1 . The synthesized beam is FWHM = 0751 x 0742 
and P.A. = -5°. (b) Map restored with a l'/5 circular beam shown in 
the bottom right corner. Thin contours: Map of the 3.5 cm continuum 
from the AC 504 dataset (Raga et al. 2000). The rms noise is 0.009 
mjy beam -1 . The contour levels are 1, 2, 4, 8, and 16 times 0.04 mjy 
beam -1 . The synthesized beam is FWHM = 0731 x 0724 and P.A. 
= -10°. Outflow knots are labeled following the convention used by 
Curiel et al. (1993). 



6.9 mm observations in D-array). While the maps from the AC 
504 dataset were more useful in identifying the sources in the 
6.9 mm maps, the flux densities from the AC 563 dataset were 
used in the analysis of source spectra. 

3. RESULTS 

Figure la shows the 6.9 mm continuum map of the SMM 
1 region from the observations in 2001. The central source 
(source 1) was clearly detected and has an extrusion toward the 
northwest. Examinations of the CLEAN components showed 
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Fig. 2. — Natural-weight maps of the 6.9 mm continuum toward the 
SMM 1 region. The 6.9 mm sources are labeled, (a) Map from the 
observations in 2004 March. The rms noise is 0.09 mJy beam" 1 . The 
contour levels are 1, 2, 4, 8, and 16 times 0.3 mjy beam -1 . Dashed 
contours are for negative levels. Shown in the bottom right corner is 
the synthesized beam: FWHM = 0769 x 0754 and P.A. = 26°. (b) Map 
from the 2001 and 2004 datasets combined. The rms noise is 0.05 mjy 
beam -1 . The contour levels are 1, 2, 4, 8, 16, and 32 times 0.15 mjy 
beam -1 . Shown in the bottom right corner is the synthesized beam: 
FWHM = 0788 x 0775 and P.A. = 15°. The straight line near the bot- 
tom corresponds to 500 AU at a distance of 260 pc. Plus signs: Position 
of the 3.5 cm sources shown in Fig. lb. 



that this extrusion is well separated from source 1. To see it 
clearly, a map was made with a restoring beam smaller than the 
usual synthesized beam (Fig. lb). While details of this "super- 
resolution" map should not be overinterpreted, it shows that the 
extrusion is an unresolved source (source 2) located close to the 
outflow axis and associated with the outflow knot E. Source 1 
seems to be elongated in the northeast-southwest direction. 

Figure 2a shows the 6.9 mm continuum map from the obser- 
vations in 2004. For source 1, the peak position of the 6.9 mm 
source agrees with that of the 3.5 cm source within 071 as well 
as with those of the H2O and OH masers present in the region 
to within 1" (Rodriguez et al. 1989; Moscadelli et al. 2006). 
The position of the millimeter source given by Hogerheijde et 
al. (1999) is ~2'/2 west of the 6.9 mm source position. The rea- 
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Table 1 

Millimeter Continuum Sources in the Serpens SMM 1 Region 



Source 


Name 


Peak Position 


Peak Flux 


Total Flux 






QJ2000.0 


<5j2000.0 


(mjy beam -1 ) 


(mjy) 


1 


SMM la 


18 29 49.80 


01 15 20.6 


7.67 ± 0.05 


13.1 ±0.4 


2 


knot E 


18 29 49.63 


01 15 22.3 


0.9 ±0.2 


unresolved 


3 


SMM lb 


18 29 49.69 


01 15 21.2 


0.47 ± 0.05 


1.10 ±0.09 



NOTES. — Units of right ascension are hours, minutes, and seconds, and units of dec- 
lination are degrees, arcminutes, and arcseconds. Flux densities at 6.9 mm are corrected 
for the primary beam response. For SMM la and lb, the flux densities are from the D- 
and C-array data combined (Fig. 26). For knot E, it is from the D-array data only (Fig. 
lb). Total flux of each source was measured in a circumscribed rectangle of 2a contour. 



son for this disagreement is not clear. We presume that source 
1 detected by us and the source detected by Hogerheijde et al. 
(1999) are actually the same object. Source 2 of the 2001 epoch 
was undetected in the map of the 2004 epoch. If its flux were 
unchanged, it should have been detected at a 10ct level, where 
a is the rms noise. Therefore, source 2 must be time variable. 
Source 3 is associated with a 3.5 cm continuum source. Though 
Curiel et al. (1993) named it knot F, it is not exactly located on 
the outflow axis. We will discuss the nature of this source in 
Section 4.2. 

Figure 2b shows the map from the D- and C-array datasets 
combined. The bright and compact central source is surrounded 
by an extended emission structure. The extended source is 
clumpy and elongated in the northeast-southwest direction with 
a full length of ~5". The flux of source 2 in Figure 2b is not 
very meaningful because it is time variable. 

Table 1 lists the parameters of the 6.9 mm continuum 
sources. In the discussion below, sources 1, 2, and 3 will be 
referred to as SMM la, knot E, and SMM lb, respectively. 

4. DISCUSSION 
4.1. SMM la 

Comparison of the centimeter and millimeter continuum 
maps clearly shows that SMM la is the driving source of the 
bipolar radio jet. To understand the origin of the 6.9 mm emis- 
sion, the spectral energy distribution (SED) of SMM la was 
constructed from flux measurements with radio interferometric 
observations (Fig. 3). The 6.2 cm, 3.5 cm, and 6.9 mm flux den- 
sities are from this work, and the 3.4/3.2 and 2.7 mm values are 
from Hogerheijde et al. (1999). The 1.4 mm flux of Hogerheijde 
et al. (1999) was not included in the plot because the minimum 
uv length of the observations was a few times larger than those 
of the other observations, which may result in resolving out con- 
tributions from large-scale structures. The uncertainties shown 
in Figure 3 include the contribution from the uncertainties of 
the absolute flux scale (2-10%). 

In the centimeter-millimeter wavelength range, the SED of 
a low-mass young stellar object usually has two components: 
dust emission from a disk/envelope system with a large (> 2) 
spectral index and free-free emission from a thermal radio jet 
with a small (< 1) spectral index (Reynolds 1986; Anglada et 
al. 1998). The SED of SMM la was fitted with a sum of two 
power-law components, F v — F\v ai ± F 2 v a2 , where F v is the 
flux density, v is the frequency, and a is the spectral index. The 
best-fit spectrum has spectral indices of ct\ = 3.82 ±0.11 and 
a 2 = -0.02 ± 0.14 and suggests that about 90% of the 6.9 mm 
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Fig. 3. — Spectral energy distribution of SMM la (filled circles), SMM 
lb (filled squares), and knot E (open circles). Solid curve: Best-fit two- 
component power-law spectrum of SMM la. Dotted straight lines: 
Each component of the fit. Filled triangle: The upper limit (3a) for 
SMM lb at 6.2 cm. 



flux comes from dust. 

With the SED of dust component, the mass of molecular gas 
can be estimated. The mass emissivity given by Beckwith & 
Sargent (1991) is assumed, 



k v = 0.1 — 



cm 2 g \ 



(1) 



where vq = 1200 GHz, and (3 is the opacity index ((3 ~ a — 2, in 
the millimeter range). If the emission is optically thin, the mass 
can be estimated by 



M 



F V D 



k v B v (T d ) 



(2) 



where D is the distance to the source, B v is the Planck function, 
and T d is the dust temperature. 
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Fig. 4. — Images of the 4.5, 5.8, 8.0, and 24 /im bands toward the SMM 1 region from the Spitzer data archive. The contour levels are in arbitrary 
linear sequences. See Harvey et al. (2007) for the flux scale. Squares: Position of the 6.9 mm sources SMM la/b. Plus signs: Position of the 3.5 cm 
sources shown in Fig. lb and VLA 8 (near the eastern edge). Long arrows: Proper motion vectors of the bipolar radio jet ((X'12 yr _1 ; Curiel et al. 
1993). Short arrow: Proper motion vector of the H2 jet (~0'.'03 yr _1 ; Hodapp 1999). 



Assuming D = 260 pc and T d = 36 K (Straizys et al. 1996; 
Larsson et al. 2000), the mass of SMM la is M = 8 ± 3 M , 
which includes the inner protostellar envelope and the disk. If 
the dust is colder (see Section 4.2), the mass can be larger. For 
comparison, the mass derived by Hogerheijde et al. (1999) is 
~4 Mq (when scaled to the same distance), and the mass of the 
envelope within an 80" beam is 15.4 Mq (Enoch et al. 2007). 
The large mass of the disk/envelope and the powerful thermal 
radio jet corroborate the suggestion that SMM la is an actively 
accreting Class protostar. 

4.2. SMM lb 

SMM lb is located l'/8 (470 AU at a distance of 260 pc) 
away from SMM la, and they were clearly separated only in the 



C-array map of the 2004 epoch. It has a steeply rising spectrum 
(a w 2.04) between 3.5 cm and 6.9 mm (Fig. 3). The SED 
suggests that a significant part of the 6.9 mm flux may be from 
dust, which would imply that SMM Ibis not an outflow knot but 
a dense core of molecular gas. If the mass is simply scaled from 
SMM la using the flux ratio as a very rough estimate, the mass 
of the molecular gas associated with SMM lb is ~0.7 Mq. 

The Spitzer images provide further useful information (Fig. 
4). While the 4.5-8.0 fim images show the outflow in the over- 
all structure (Harvey et al. 2007), the brightest mid-IR source 
(source 141) coincides with SMM lb, and SMM la was un- 
detected. The peak position of the 24 /im source is situated in 
the middle between SMM la and lb. In other words, SMM lb 
dominates the mid-IR part of the SED, SMM la dominates in 
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the millimeter range, and they may be comparable in the far-IR 
flux. 

As shown in Figure 16, SMM lb is comparable to other 
clumpy structures around SMM la (e.g., the clump l'/5 north- 
east of SMM la) in terms of intensity. However, SMM lb is 
different from the other clumps in several respects. First, only 
SMM lb was clearly detected in the C-array map (Fig. la), 
which suggests that SMM lb has a compact structure. Second, 
SMM lb has counterparts in other wavelengths: knot F in the 
centimeter continuum images and source 141 in the IR images. 
Third, the other clumps seem to belong to an extended struc- 
ture elongated in the direction perpendicular to the bipolar jet, 
which suggests that it is a disk-like fiat structure around SMM 
la, while SMM lb is located away from this extended structure. 
The most likely interpretation is that SMM lb is yet another 
young stellar object in this region. In other words, SMM 1 may 
be a protobinary system. 

The SMM la/b system illustrates the importance of high- 
resolution imaging. The SED of SMM 1 in previous studies 
(e.g., Larsson et al. 2000) is actually a combination of the SEDs 
of the two sources. While an accurate SED of each source is not 
available unless the system can be resolved in the far-IR band, 
the dust temperature of SMM la may be lower than what was 
derived previously. SMM lb is probably less deeply embedded 
than SMM la but still very young as it was not clearly detected 
in the 3.6 /im band and undetected in shorter wavelengths. 

In the classification scheme of multiple young stellar ob- 
jects using the binary separation, as suggested by Looney et al. 
(2000), SMM 1 belongs to the common-envelope type. It will 
be interesting to image the system with a higher angular reso- 
lution to see the relation between the two binary components, 
such as the alignment of the disks. 

It is interesting to note that there is an H 2 jet that seems to 
be driven by SMM 1 (Hodapp 1999). It can also be seen in the 
IR images (Fig. 4). Since the direction of the H2 jet has a posi- 
tion angle different from that of the bipolar radio jet by ~25°, 
the two jets may have separate driving sources. SMM lb is a 
probable driving source of the H 2 jet. 

4.3. KnotE 

Knot E has a flat (a « 0.09) spectrum (Fig. 3), which sug- 
gests that the 6.9 mm flux is mostly free-free emission from the 
thermal jet. The nondetection at the 2004 epoch suggests that 
knot E is fading rapidly. Such a fading has been observed in 
other radio jets such as the HH 80-81/GGD 27 jet (Martf et al. 
1998). 

4.4. VIA 8 

VLA 8 is a centimeter continuum source located 23" north- 
east of SMM 1 (Eiroa et al. 2005). It was detected in the 6.2 and 
3.5 cm images as an unresolved source with flux densities of 
0.481 ± 0.010 and 0.442 ± 0.011 mJy, respectively. The peak 
position is a 20 oo = 18 h 29 m 51!04, <5 2 ooo = 01°15'33'/9. VLA 8 
was undetected at 6.9 mm (< 0.24 mJy beam -1 ). The SED sug- 
gests that the centimeter continuum flux comes from free-free 
emission and that there is no substantial emission from dust. 
The nature of VLA 8 is not clear. 
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